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^ : ABSTRACT 

in 

Q ■ The stellar core formation and high speed jets driven by the formed core 

are studied by using three-dimensional resistive magnetohydrodynamical (MHD) 
^ nested grid simulations. Starting with a Bonnor-Ebert isothermal cloud rotating 

Qh' in a uniform magnetic field, we calculate the cloud evolution from the molecular 

O ' cloud core {uc = 10^ cm~^, rc = 4.6 x 10^ AU) to the stellar core {uc — 10^^ cm~^, 

' rc — I Rq), where Uc and rc denote the central density and radius of the objects, 

■ respectively. We resolve cloud structure over 7 orders of magnitude in spatial 
extent and over 17 orders of magnitude in density contrast. For comparison, we 

^ ' calculate two models: resistive and ideal MHD models. Both models have the 

■ same initial condition, but the former includes dissipation process of magnetic 
field while the latter does not. The magnetic fluxes in resistive MHD model 
are extracted from the first core during 10^^ cm^'^ < Uc < 10^^ cm~'^ by Ohmic 
dissipation. Magnetic flux density of the formed stellar core [uc — 10^° cm~'^) in 
resistive MHD model is two orders of magnitude smaller than that in ideal MHD 
model. Since magnetic braking is less effective in resistive MHD model, rapidly 
rotating stellar core (the second core) is formed. After stellar core formation, 
the magnetic field of the core is largely amplified both by magneto-rotational 
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instability and the shearing motion between the stellar core and ambient medium. 
As a consequence, high speed (~ 45 kms~^) jets are driven by the second core, 
which results in strong mass ejection. A cocoon-like structure around the second 
core also forms with clear bow shocks. 

Subject headings: ISMrjets and outflows — ISM:magnetic flelds — MHD — stars: formation 

1. Introduction 

Outflows and jets are ubiquitous in star-forming regions and closely related to the ac- 
cretion process, angular momentum transportation and magnetic field in the cloud. Recent 
observations indicate these flows are composed of multiple components in space and velocity 
(Pyo et al. 2003, 2005). Molecular outflows observed by CO emission have speed of ~ 10 
kms~^ and wide opening angle (e.g., Belloche et al. 2002), while jets observed by SiO 
and Fell emissions have speed of 100 kms~^ and well coUimated structure (e.g., Hirano et 
al. 2006). These observations suggest that different flows which have different sizes, flow 
speeds, and opening angles are driven by similar mechanism but with different conditions, or 
at different evolutional phases. Although internal structure and rotational motion inside the 
outflows and jets arc observed by recent high-resolution telescopes (Bacciotti et al. 2002; 
Choi 2005), driven points of these flows arc not observed because these regions arc embedded 
in the dense cloud core and located in the very vicinity of the protostar where we can not 
spatially resolve. Thus, the driving mechanism of outflows and jets are not understood well. 
Since the effect of radiation pressure is small in low mass stars, these flows are supposed 
to be driven by the Lorentz force. The magnetic energies observed by Zeeman effect mea- 
surements are comparable or larger than the thermal energy (Crutcher 1999), implying that 
magnetic effects are important for the evolution of molecular clouds. Thus, numerical MHD 
simulations are needed in order to understand the formation processes of protostars, jets and 
outflows. 

Assuming spherical symmetry, the evolutions from the molecular cloud to the stellar core 
are calculated by many authors (e.g., Larson 1969; Winkler & Newman 1980a,b; Masunaga 
et al. 1998; Masunaga & Inutsuka 2000), and they found density, velocity and temperature 
structures of the cloud and stellar core. Bate (1998) and Whitchousc & Bate (2006) have 
calculated stellar core formation from the molecular cloud core in their three dimensional 
SPH simulations. In their calculations, the spiral or bar structures are appeared at high 
density, because the bar mode instability is induced in rapidly rotating cores (Durisen et al. 
1986). Since angular momentum is effectively removed due to non-axisymmetric structure, 
the central region contracts and stellar core forms. However, angular momentum is largely 
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removed by the magnetic braking (Basu & Mouschovias 1994; Machida et al. 2005a) and 
by driving outflows (Matsumoto & Tomisaka 2004; Machida et al 2004, 2005b). Therefore, 
they overestimate the angular momentum of the formed core, because they ignore the mag- 
netic field in the cloud. The evolution of magnetized clouds are investigated by Tomisaka 
(1998, 2000, 2002) in his two dimensional nested grid simulation. In his calculation, low- 
velocity (outflow; V ~ 2kms^^) and high-velocity (jet; v ~ 30 kms~^) flows are driven 
by first core and second core, respectively However, he adopt ideal MHD approximation, 
which is valid in low-density gas region {n ^ lO^^cm"^), however, not valid in high-density 
gas region (n ^ 10^^ cm~^). Nakano et al. (2002) found significant magnetic fiux loss occurs 
during lO^^cm"^ ^ n ^ lO^^cm"^ by Ohmic dissipation. Therefore, Tomisaka (1998, 2000, 
2002) overestimate the magnetic fiux of the cloud especially in high-density gas region. 

In this paper, we calculate cloud evolution from cloud core {ric — 10^ cm""^, Vc — 4.6 x 10^ 
AU) to stellar core phase {ric ^ 10^^ cm~^, Tc ~ 1 Rq) using three dimensional resistive MHD 
nested grid method and study the formation process of high-speed jets. We also calculate 
cloud evolution using ideal MHD approximation, and compare the resistive MHD model with 
the ideal MHD model. 



2. Model and Numerical Method 

Our initial settings are almost the same as those of Machida et al. (2006). To study 
the cloTid cvohition, wc use the three-dimensional resistive MHD nested grid code. We solve 
the resistive MHD equations including the self-gravity: 

^ + V-(pt;) = 0, (1) 

p-^ + p{v ■ V)v = -VP - — S X (V X B) - pV0, (2) 
dB 

— ^Vx{vxB) + TjV^B, (3) 

VV = 47rG'p, (4) 

where p, v, P, B, r] and (p denote the density, velocity, pressure, magnetic fiux density, 
resistivity and gravitational potential, respectively. To mimic the temperature evolution 
calculated by Masunaga & Inutsuka (2000), we adopt the piece- wise polytropic equation of 
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state. 



cIp P < Pc, 

Cipc [fj Pc<P< Pd: 



where = 190ms~\ pc = 1.92xl0~^^g cm'^ {uc = 5 x 10^° cm-^), pd = 3.84 x 10~^ g cm~3 
(ud = lO^^cm"^), and pe = 1.92 x 10~^g cm^^ (rie = 5 x 10^'^ cm"^). In this paper, a 
spherical cloud with critical Bonnor-Ebert (Ebert 1955; Bonnor 1956) density profile having 
Pc,o = 3.84 X 10~"'^^g cm~^ {nc,o = 10^ cm^^) of the central (number) density is used as the 
initial condition. Initially the cloud rotates rigidly {Qq = 1-02 x 10^^'^s^^) around the 2;-axis 
and has uniform magnetic field {Bq = 212/iG) parallel to the 2;-axis (or rotation axis). To 
promote the contraction, we increase the density as by 30 % from the critical Bonnor-Ebert 
sphere. 

We adopt the nested grid method (for detail, see Machida et al. 2005a) to obtain high 
spatial resolution near the center. Each level of rectangular grid has the same number of 
cells (= 64 X 64 X 32), but cell width h{l) depends on the grid level I. The cell width is 
reduced to 1/2 with increasing grid level {I — > / + 1). The highest level of grids changes 
dynamically. We start the calculation with 4 grid levels {I — 1,2,3,4). Box size of the 
initial finest grid / = 4 is chosen equal to 2Rc, where Rc denotes the radius of the critical 
Bonnor-Ebert sphere. The coarsest grid (/ = 1), therefore, has a box size equal to 2* R^.. A 
boundary condition is imposed at r = 2*^ i?c, where the magnetic field and ambient gas rotate 
at an angular velocity of Oq (for detail see Matsumoto & Tomisaka 2004). A new finer grid 
is generated whenever the minimum local Jeans length Aj becomes smaller than 8/i(/max), 
where h is the cell width. The maximum level of grids is restricted to /max = 21. Since the 
density is highest in the finest grid, generation of new grid ensures the Jeans condition of 
Truelove et al. (1997) with a margin of a safety factor of 2. 

We calculate two models: (a) resistive and (b) ideal MHD models. Both models have the 
same initial condition shown above. The former include resistive term {rjV'^B ) in induction 
equation (eq. [3] ), while the latter does not. 
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3. Results 

The molecular gas obeys the isothermal equation of state with temperature of ~ 10 
K until nc ~ 5 X 10^°cm~^ ( isothermal phase), then cloud collapses almost adiabatically 
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(5 X 10^° ^ 10^^; adiabatic phase) and quasi-static core (i.e., first core) forms during the 
adiabatic phase (Larson 1969; Masunaga et al. 1998). In our calculations, the first core forms 
both in resistive and ideal MHD models when the central density reaches Uc — 8x 10^^ cm~^ 
and these cores have radius of ~ 10.6 AU. The magnetic flux is expected to be removed from 
the first core during the adiabatic phase (Nakano et al. 2002). We estimate the resistivity 
(?7) according to Nakano et al. (2002) and represented by a thick line in upper panel of 
Figure 1 which has a peak at n ~ lO^^cm"^. We also estimate the magnetic Reynolds 
number {R = vj Xj r)'^; dotted line in upper panel of Fig. 1) using the free-fall velocity (vf) 
and Jeans length (Aj). The dotted hne indicate that magnetic dissipation is effective during 
10^^ cm~^ ^ric ^ 10^^ cm~^ which corresponds to the results of Nakano et al. (2002). 

The ratio of magnetic to thermal energy within n < 0.1 and the angular velocity 
normalized by the free-fall timescale at the center are plotted against the central density in 
lower panel of Figure 1. The magnetic energy and angular velocity are almost the same during 
10^ cm~^ < ric <^ 10^^ cm~^ between resistive and ideal MHD models. The magnetic energy 
decreases in early isothermal phase (10® cm^"^ < <^ lO^cm"^), because the cloud collapses 
along magnetic field line (Machida et al. 2005a), then the magnetic energy is comparable 
to the thermal energy during 10^ cm^^ <^ lO^^cm"^. The ratio of magnetic to thermal 
energy in ideal MHD model keeps a constant value during 10® cm~'^ ^ nc ^ 10^° cm~^ and 
the magnetic field have equivalent energy to the thermal energy until final stage (the ratio 
of the magnetic to thermal energy is ~ 0.25). The angular velocities both in resistive and 
ideal MHD models decrease during isothermal phase, because cloud collapses vertically and 
magnetic braking is effective in this phase as shown in Machida et al. (2005a). The angular 
velocity in ideal MHD model begins to increase at n ~ 10^°cm~^, then it keeps an almost 
constant value (culAnGp]^^^'^ ^ 0.02) until the final stage, while it continues to increase after 
the central density reaches ric — 10^^ cm^"^ in resistive MHD model. The evolutional tracks 
of magnetic energy and angular velocity begin to depart from each other at n ~ 10^^ cm~^. 

The magnetic energy in resistive MHD model begins to decrease at Uc — 10^^ cm^^, 
and becomes about 10~^ times smaller than thermal energy in resistive MHD model at 
ric — 10^^ because magnetic field lines arc extracted from the high density region by the 
Ohmic dissipation. Conversely, the angular velocity in resistive MHD model becomes 10 
times larger than that in ideal MHD model at ric — 10^°cm~^. Since the magnetic braking 
in resistive MHD model is less effective than that in ideal MHD model for weak magnetic 
field, the angular momentum in resistive MHD model is removed slightly compared with ideal 
MHD model. After central density reaches ric — lO^^cm"^, the equation of state becomes 
soft (see equation [5]) refiecting the dissociation of hydrogen molecules at T ~ 2 x 10^ K, and 
collapses rapidly, i.e., the second collapse begins. By this epoch, the central temperature 
becomes so high that the thermal ionization of Alkali metals reduces the resistivity and so 
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that Ohmic dissipation becomes ineffective. Thus, the magnetic field becomes strong again 
as central region collapses. The second core (cf. Larson 1969) formed at ric — 10^^ cm~^ 
both in resistive and ideal MHD models. The magnetic field strength increases rapidly after 
the second core formation epoch (n ^ 10^^ cm^"^), because the shearing motion between the 
second core and ambient medium amplifies the toroidal magnetic field around the second 
core. 

Figure 2 shows the time sequence of the cloud in resistive MHD model after the second 
core formation (n^ ^ 10^^ cm^'^). Panels in this figure show only grids of 1=17 and 18 levels 
to resolve jet-like structure, while we have resolved the deep interior at the central region 
with the grids of / = 1 — 21 levels. A thin disk is formed inside the first core when the 
central density reaches ric ^ 1.2 x 10^^ cm~^, and then a spherical core is formed within the 
disk. The disk has 19.4 Rq of radial scale and the spherical core has 5.5 Rq of radius. The 
mass of disk and core, respectively, increase up to 1.6 x 10~^ Mq and 7.2 x 10~^ Mq at the 
end of the calculation. We stopped calculation at 10.86 days after the second core formation 
(ric = lO^^cm"^), because the time step becomes extremely small (At ~ 25 second !) for 
the Alfven speed being high in the region just outside of the disk. The green contour in 
each panel denotes the iso- velocity curves representing the jet region in which the gas is 
outfiowing from the core. The gas is outfiowing from the spherical core with wide opening 
angle (~ 45°) and 25 km s~^ of the maximum speed in Figure 2a. The jet expands to 
the vertical direction and cocoon- like structure is formed as shown in Figure 2b. At the 
calculation end, 7.2 x 10~^ Mq of the gas is outflowing with 41 km s~^ of the maximum 
speed. The average accretion and outflow rate which are calculated after jet appeared are 
3.4 x 1Q~^ Mq yr~^, and 2.2 x 1Q~^ Mq yr~^, respectively. The bow shock structures are 
clearly seen at 2; = ±0.24 AU in Figure 2c. 

Figure 3 shows the cloud structure and magnetic fleld lines in the bird's-eye view at the 
same epoch as Figure 2c. The structure of the jet is shown by blue iso- velocity surface. The 
jet is coiled by the magnetic fleld (stream lines). The red iso-surfaces above and below the 
second core inside blue iso-surface means the high-density region (ric > 10^® cm^"^) which 
indicates the strong mass ejection from the center that occurred in the past. The ejected 
mass is also coiled by the magnetic field. 

4. Discussion 

We start with a slowly rotating and strongly magnetized cloud in this paper. The ratio 
of rotational to gravitational energy is 4 x 10~^, while the magnetic energy is equivalent 
to the thermal energy in the initial cloud. The outfiow driven by the first core is seen in 
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Tomisaka (2002), Matsumoto & Tomisaka (2004), and Machida et al. (2004, 2005b), while 
these outflows does not appear in this model. This is because considerably slowly rotating 
cloud is assumed initially. However, we confirmed that the outflow is driven also by the first 
core in a more rapidly rotating cloud even if resistive MHD is adopted. 

The high-speed jet (~ 40 kms~^) is driven by the second core both in resistive and ideal 
MHD model. We have also calculated other models with various sets of parameters using 
both of resistive and ideal MHD methods. The jet always appears in ideal MHD models, 
while jet does not appear in some resistive MHD models by ~ 10 days after the second core 
formation. As shown in §3, the magnetic energy in resistive MHD model is much smaller 
than the thermal energy when the second core is formed. If the magnetic energy is larger 
or comparable to the thermal energy, jets are easily driven from the formed core by the 
magneto-centrifugally wind mechanism (Blandford&Payne 1982). However, this mechanism 
is inefficient when the magnetic energy around the protostar is extremely small. In our 
resistive MHD model, the second core has a small magnetic energy. Instead, the magnetic 
field in the formed disk around the central core is gradually amplified. We can see from the 
density structure in the z = plane projected into bottom wall of Figure 3 that the spiral 
patterns are formed. When the central density exceeds ric > 10^^ cm^^, non-axisymmetric 
(or spiral) patterns begin to grow inside the disk. It is considered that these patterns are 
caused by 'bar model instability' (e.g., Durisen et al. 1986) or magneto- rotational-instabihty 
(Balbus & Hawley 1991). Because of these instabihties, the second core accretes not only 
gas but also magnetic flux from the surrounding disk. This results in the increase of the 
magnetic pressure that is sufficient to launch the jet. The high speed jet in resistive MHD 
models seems to require these instabilities. Note that we had to stop the calculations in all 
models at ~ 20 days after the second core formation due to the extremely small time step. 
Thus, it remains to determine whether or not and when the jet is driven ultimately in the 
models that have not shown it so far. However, jets may appear in models in which jet does 
not appear if we calculate further. The detailed analysis of the mechanism and conditions 
for driving the jets will be investigated in our subsequent papers. 

We have greatly benefited from discussion with Prof. K. Tomisaka, Dr. K. Saigo and 
Dr. Y. Kato.We also thank T. Hanawa for contribution to the nested grid code. Numerical 
calculations were carried out with a Fujitsu VPP5000 at the Astronomical Data Analysis 
Center, the National Astronomical Observatory of Japan. This work was supported partially 
by the Grants-in-Aid from MEXT (15740118, 16077202, 16740115). 
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Fig. 1. — Upper panel: Resistivity [r]; left axis) and magnetic reynolds number (R; right 
axis) as a function of the central number density. The magnetic field well couple with the gas 
in "coupled" region, while magnetic field decuples from the gas in "decoupled" region. Lower 
panel: The magnetic (B'^/Stt) to thermal energy (c^p) ratio (left axis ) within n < O.lric and 
angular velocity (i7) normalized by free-fall timescale {[AnGpcY^'^) center (right axis) 

against central number density. The gas begins to ionized when number density reaches 
n ~ 10^^ cm~^. 
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(=2.3458948x10= [Myr] n„„= 1.6x10"' [cm '] (=2.3458957x10= [Myr] iw= i 7x10^' [cm"'] (=2.3458960x10= [Myr] n^,,= 1 .8x10'' [cm-^] 
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Fig. 2. — The density (color-scale and white contours) and velocity distribution (arrows) on 
the cross-section in the y = plane. Panels (a)-(c) are snapshots at different stage, but the 
same grid level (/ = 17 and 18). Green thick lines are velocity contours of outflow {vout > 3 
km s~^; thick line; Vout >20 km s~^; dashed line). 



Fig. 3. — Bird's-eye view of Z = 17 grids in Fig. 2 c. The structure of high-density region 
(n = lO^^cm"^; red iso-density surface), velocity vectors of outflow (thick arrows), and 
magnetic field lines are plotted. The structure of the jet is shown by blue iso-velocity surface 
in which the gas is outflowing from the center. The density contours (false color and contour 
lines), velocity vectors (thin arrows) are projected in each wall surface. 



